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Abstract. The aim of this paper is twofold. First of all, we study the behaviour 
of the lowest eigenvalues of the quantum anharmonic oscillator under influence of an 
external field. We try to understand this behaviour using perturbation theory and 
compare the results with numerical calculations. This brings us to the second aim of 
improving the method used to carry out the numerical calculations. 



1. Introduction 

The quantum anharmonic oscillator is an important model in physics. Unfortunately 
it can't be solved in closed form. In the past, intensive efforts were carried out to 
develop reliable methods to solve the model numerically. Many papers study only the 
double well potential [H El [3] . Other papers deal with all sorts of quantum anharmonic 
oscillators (see, e.g. [1] and references therein). The main point of these papers is to 
develop numerical methods to solve the eigenvalue problem of anharmonic oscillators as 
accurately as possible. Selecting the best method out of all these different proposals is 
a hard task. 

The problem of the quantum anharmonic oscillator in an external field has not 
yet been studied in detail, although it has interesting applications. In [5] one studies 
numerically chains of coupled quantum anharmonic oscillators. Such chains are used in 
condensed matter as one dimensional models of crystals and are a first interesting step 
towards quantum field theory. In this context, coupled anharmonic oscillators are used 
to study displacive and order-disorder phase transitions (see e.g. [6]). The anharmonic 
oscillator in an external field appears in the mean field treatment of a chain of coupled 
anharmonic oscillators when the interactions between the oscillators are replaced by an 
effective external field. 

In Section 2, we will follow the lines of ffj to calculate the eigenvalues of the 
anharmonic oscillator with a general polynomial potential. In Section 3, we apply the 
general results obtained in Section 2 to the anharmonic oscillator in an external field. 
In Section 4 and 5 we study the lowest eigenvalues of respectively a shallow and a deep 
well. The final section gives a short discussion of the results. 
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2. Method 

We use a method which is based on two papers [21 E]- In [7J, Jafarpour and Afshar 
calculate the energy eigenvalues of the anharmonic oscillator with a general polynomial 
potential. To obtain these eigenvalues they define a 'normalised squeezed vacuum state' 
and use it to build up an orthonormal basis. With a variational approach, this basis is 
optimised. Then it is used to calculate the matrix representation of the Hamiltonian. 
Finally this matrix is diagonalised numerically to obtain the lowest energy eigenvalues. 
Following the lines of |7j we calculate the eigenvalues of the general Hamiltonian 

H=^p'+j:w- (1) 

However, we adapt the method of [3| at two important places. First, we use the well 
known basis of the harmonic oscillator in stead of the transformed basis introduced in 
[7]. This makes the method more transparent. Second, we use a more general method 
for the optimisation of the choice of basis. 

1) Our modification of the method starts from the observation that the 'normalised 
squeezed vacuum state' of [7| is just the ground state of a harmonic oscillator 
with suitably chosen parameters. Indeed, the transformation from one harmonic 
oscillator to another one, with different mass m and frequency uo, is a special case of 
a Bogoliubov transformation. The generic case of the Bogoliubov transformation 
is not used in |7]. The basis functions used by the authors are the well known 
eigenf unctions of the harmonic oscillator, be it with optimised value of the 
parameters m and c^o- The method of [7] can therefore be simplified by omitting 
the use of Bogoliubov transformations. 

2) In the literature, other methods [H [2] are found than that of [7], to optimise the 
choice of orthonormal basis. We will discuss some of these methods and use the 
one best suited for the problem at hand. 

A more detailed description of the method now follows. 

The operators a and are the annihilation and creation operators of the ordinary 
harmonic oscillator with mass m and frequency Uq. With these operators and the 
notation muJor"^ = h Hamiltonian ([T]) becomes 



H 



(a -at) +EA.(^-^j (a + at)\ (2) 

An ordered expression for (a + a^) is derived in |7j. Using this expression the general 
Hamiltonian can be written as 

[a;J means rounding the value of x to the lower integer. 



The quantum double well anharmonic oscillator in an external field 



3 



The method proceeds in three steps. First, the expectation value of the Hamiltonian 
([3]) is calculated in some eigenstate of the harmonic oscillator (the choice of eigenstate 
will be discussed later on). Second, the expectation value is minimised by varying the 
parameter UQm, or equivalently r, of the harmonic oscillator. Third, the eigenfunctions 
of the harmonic oscillator (now with fixed parameter) are used to obtain the matrix 
representation of the Hamiltonian H. 

The expectation value in an arbitrary eigenstate, labeled with parameter t, is given 

by 

H \A) = ^i2t + 1) + E 2V2+.fe,((,/2_'^),).(,_,/2 + fc)! - 

After minimising this expectation value {d {iptl H \tpt) /i9r^=0), the following equation 
is obtained 

^ J. z\,rl^i\t\ 

°~ 4^^^^ + 1) + ^X.^ 2^n+k+ik\{{i/2 - k)\f{t - i/2 + A;)! ' ^' 

This equation gives one condition to calculate the value of two parameters (ro and t). In 
practice we will chose a value for the parameter t, and solve the equation to obtain the 
value of ro, the optimal value of r. In general this equation has to be solved numerically. 
In the following section, we return to the problem of choosing the parameter t. 

The eigenfunctions of the harmonic oscillator with r = r^ are now used as a basis. 
The matrix elements of Hamiltonian ([3]) become 



AmrQ 



{2t + l)6s,t - ^t{t - l)5s,t-2 - ^{t+l){t + 2)5s,t+2 



L*/2J i-2k 



^ § S S ^^^^^j'(t-j)!(^-2A:-jV^'*^^-^-^^- 

Note that the last part equals zero, if j > t and i — 2k — j > s. The matrix is infinite 
dimensional. Truncating the basis to N basis functions results in an N dimensional 
matrix. After diagonalising this dimensional matrix, the low-lying eigenvalues of the 
Hamiltonian are found. 



t\{t + i - 2j ~ 2k)\ 



3. Anharmonic oscillator in an external field 

From now on we restrict ourselves to the quantum anharmonic oscillator in an external 
field. The Hamiltonian is given by 

with a < and /? > 0. The potential V (Q) is a double well potential. The extra term 
—pQ favours one of the two wells. 
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Rewriting this Hamiltonian to the form of ([T]) gives / = 4, Aq = A3 = 0, Ai = —p, 
A2 = a/2 and A4 = (3/4. Then equation ([5]) has to be solved, to obtain the value of rg. 
In this case, the equation reduces to 



This is a cubic equation in Tq which has only one real and positive solution (remember 
a < 0). 

For the double well potential different choices of the parameter t are proposed in 
the literature. None of these choices has a physical interpretation. The only reason to 
prefer one choice above another is a faster convergence of the numerical results as a 
function of A^. 

The condition derived in [7] corresponds with the choice t = 0. This is not the 
optimal choice. The authors are aware of this and do not always use the calculated value 
of To to obtain the numerical data. In [1] Balsa et al. minimise the expectation value 
dl]) with respect to the parameters r and t. This way, the authors get two conditions 
for two parameters and no arbitrary choice is necessary. Problem is, as outlined in [2] 
by Bishop et al., that the calculated values of the parameters are far from the optimal 
values. 

Bishop et al. also noticed that the optimal choice of ro depends of A^. This motivates 
the choice t = N/2. This is the choice we will use to obtain our numerical results. As 
mentioned earlier, there is no physical justification for this choice. The only motivation 
is that we get a fast convergence of the numerical results as a function of A^ for every 
example we studied. 

We are interested in the behaviour of the eigenvalues under influence of an external 
field. So we have to repeat the diagonalisation of the Hamiltonian several times for 
different values of the parameter p. Note that we don't have to repeat the calculation 
of ro because equation ([8]) doesn't depend on p. 

4. Shallow well 

For the first example we choose a = —2 and (3 = 1. 

Figure (PA) shows the eigenvalue of the ground state and the first excited state 
as a function of the parameter p. If p = 0, the potential is the symmetric double well 
potential. If an extra term —pQ is added, the energy decreases, because the system 
can lower it's energy by partly choosing one of the wells. However, the energy of the 
first excited state increases for small values of p. To understand this, we look for the 
response of the wavefunctions to the parameter p. 

Figure ([IB), shows the wavefunction of the ground state and the first excited state 
in the position representation, with p = 0. Clearly, the moduli of the wavefunctions 
of the ground state and of the first excited state are both symmetric. If p 7^ this 
symmetry is destroyed. Figure ([2JA), shows the wavefunctions of the ground state and 
the first excited state in the position representation, with p = 0.2. We see that the 
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Figure 1. (A) First two eigenvalues of the anharmonic oscillator in an external field 
as a function of the parameter p, with N =^ 30. (B) Position representation of the 

wavefunctions of the ground state ( ) and the first excited state ( ) of the 

anharmonic oscillator without external field {p = 0), with N = 30. 




Figure 2. (A) Position representation of the wavefunctions of the ground state ( ) 

and the first excited state ( ) of the anharmonic oscillator in an external field 

{p = 0.2), with N = 30. (B) Comparison of the result of equation (|10p ( ), 

equation (|12p (• ) and the original numerical data ( ) for the ground state. We 

used following values ci — —3.3921, a ~ —0.8531, ui — ci/a — 3.9762 and N = 40. 



wavefunction of the ground state has a large peak for positive values of q. Because of 
the orthogonality of the wavefunctions, the wavefunction of the first excited state must 
have it's largest contribution for negative values of q. As a consequence, the energy of 
the ground state will decrease and the energy of the first excited state will increase. 

Let us now focus on the response of the ground state eigenvalue to an external field. 
We consider the external field as a perturbation of the anharmonic oscillator. According 
to non-degenerate perturbation theory to second order, the energy of the ground state 
is given by 



m(^0) 



Eo(0) - E^O) 
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Table 1. The values of Tq are calculated with equation ^ and the choice t — N/2. 
In this table, the values of Tq are rounded but for the numerical calculations we used 
the exact solution of this equation. (A) The results of a numerical calculation of 
i?o(0), Ei(0) and ci, with p = and for different dimensions of the Hamiltonian 
matrix. (B) Fitting parameter for different dimensions of the Hamiltonian matrix. 







(A) 






(B) 




N 




Eo 


El 


Cl 


a 




10 


0.59 


-0.299479413549 


0.046558837188 


-3.390313017166 


-0.930574586052 


20 


0.45 


-0.299521364979 


0.046371082733 


-3.392128162181 


-0.930^ 


^91427590 


30 


0.38 


-0.299521367416 


0.046371082228 


-3.392128193573 


-0.930^ 


^91443755 


40 


0.34 


-0.299521367416 


0.046371082228 


-3.392128193573 


-0.930^ 


^91444167 



The modulus of \ipo{0)) is symmetric in q. As a consequence, the hnear term in 
expression (Q vanishes and the leading order terms of expression (Q are 



Eoip) - ^oW+J^- '^-^^Q^l'gy;' = EoiO) + cip\ (10) 

because i?o(0) — -Ei(O) << Eq{^) — Ei{Q) for z > 1. The results of a numerical calculation 
of Eq{Q), -E'i(O) and ci for different values of are given in table dHA.). As is clearly 
seen in figure ([2]B), equation (fTO!) is only valid for very small values of the parameter p. 

For intermediate values of the parameter p, the tunnel splitting between the ground 
state and the first excited state can be ignored. As a consequence, the ground state 
becomes twofold degenerate. Degenerate perturbation theory gives 

Eo{p) ^ i5o(0) -p|(^o(0)|Q |^i(0))| = Eo(0) -Pigoil • (11) 

A numerical calculation of this matrix element gives following approximate value 
igoil ~ 0.853104. 

We therefore expect that the ground state can be fitted for small and intermediate 
values of the parameter p by following equation 

Eo{p) ^ Eo(0) + aptanh(cjp). (12) 

For the limit p ^ the desired quadratic behaviour Eq{p) ^ -Eo(O) + aujp^ is obtained. 
For large p, tanh(a;p) goes to one exponentially fast. There remains a linear equation 
£;o(p) ~^o(0) + ap. 

We fix the value of uj by the condition auj = ci. The only remaining fitting- 
parameter is then a. The values of a for different dimensions of the Hamiltonian matrix 
are shown in table dlB)- The value of a ^ —0.931 is approximately the same as the value 
of — IQoil mentioned above. This means that it is a good approximation to neglect the 
tunnel splitting for intermediate values of p. 

Figure (I2]B) shows the result of equation ( ITOl) and equation ( |T2i) together with the 
original numerical data (with A^ = 40). The figure demonstrates that the non-degenerate 
perturbation theory holds for very small values of the parameter p. Equation ( JT2l) holds 



The quantum double well anharmonic oscillator in an external field 



7 




Figure 3. (A) First five eigenvalues of the anharmonic oscillator in an external field 
as a function of the parameter p. (B) Comparison of the result of equation and 

equation IjlSp ( ) and the original numerical data ( ) in the surroundings of 

pi. We used following values Qn = -0.05912, ca = -39,30905 and Q22 = -0.05902. 



for small and intermediate values of p. However, it is clear that for large values of the 
parameter p higher order terms become important. We conclude that equation (fT2|) in 
combination with the results of the perturbation theory (a = — |Qoi| and au = ci) is a 
good approximation for small and intermediate values of p. No fitting-parameters are 
necessary. 

For very large values of the parameter p, the depth of the well is asymptotically 
given by —0A725p'^^'^/X\^^. We therefore expect for very large p-values the following 
behaviour of the ground state energy 

Eoip)^A + Bp'/', (13) 

with A and B constants. This relation is in agreement with numerical calculations. 



5. Deep well 

For this second example we choose a = —4 and (3 = 1. All numerical calculations for 
this example are carried out with N = 50. 

Figure (I3]A.) shows the first five eigenvalues as a function of the parameter p. This 
figure demonstrates the repulsion of nearly degenerate energy levels. We can understand 
this from perturbation theory. Let us focus on the repulsion of the first and second 
excited state. The value of p where the energy difference between these two states is 
minimal, is denoted pi and approximately equals 0.70724. For Ei, the leading order 
terms of perturbation theory around pi are 



Eiipi + Ap) ^ Ei{pi) + Ap {Mpi)\Q \Mpi)) + i^P) 



2\{Mpi)\Q\Mpi))f 



E,{p,) - E2ipi) 

E,{p,)+QnAp + C2{Apy. (14) 
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The linear term does not vanish because the modulus of {ipn{p)\ is not symmetric for 
p ^ 0. A similar expression can be derived for E2 

E^ipi + Ap) ^ E2{pi) + Q22AP - C2(Ap)2. (15) 

The values of Ei (p) and E2 (p) are almost equal in pi. As a consequence the linear terms 
of equation (fT4l) and equation (ITS!) are also almost equal. Because of the opposite signs 
of the quadratic terms in these equations, the perturbation will lower one energy level 
and raise the other. 

A numerical calculation gives the following approximate values Qn ~ —0.05912, 
C2 ^ -39.30905 and Q22 ~ -0.05902. Figure §B) shows the result of equation ([HD 
and equation (fTSl) together with the original numerical data in the surroundings of pi. 
It's clear that the perturbation series give good approximations for the energy levels in 
the surroundings of pi. 

6. Discussion 

We study the behaviour of the lowest eigenvalues of the quantum anharmonic oscillator 
under influence of an external field. The method used to obtain our numerical data is 
based on two papers [21 [7j . We follow the lines of [7] to calculate the eigenvalues of the 
general Hamiltonian ([1]). We use the eigenfunctions of the harmonic oscillator as a basis. 
We optimise this choice of basis by varying a parameter of the harmonic oscillator, using 
a criterion inspired by [2]. Finally this matrix is diagonalised numerically. 

We apply the method to study the anharmonic oscillator in an external field. In 
[Tj, Jafarpour and Afshar considered the anharmonic oscillator as a test example. We 
were able to reproduce their results but with a quicker convergence. 

We are especially interested in the response of the energy of the ground state to 
an external field. The anharmonic oscillator in an external field is also used in [3] as 
a test example for still an other technique to calculate the eigenvalues of anharmonic 
oscillators. The authors looked also to the response of the energy levels to an external 
field, but only for intermediate values of the external field. In our opinion, there are 
three different regimes. In the first regime (small external field) , we use non-degenerate 
perturbation theory to explain the quadratic response of the energy of the ground 
state. In the second regime (intermediate external field), we use degenerate perturbation 
theory to explain the linear response of the energy of the ground state (just like in [3]). 
In the third regime (large external field), higher order terms become important. We 
propose formula (fT2l) . which is able to produce the desired behaviour in the combined 
regime of small and intermediate external fields. This formula is supported by numerical 
calculations, as is clearly visible in figure (^B)- 

In [3] , the authors also noticed the repulsion of the nearly degenerate energy levels. 
Again they only looked at intermediate external fields. They explained with degenerate 
perturbation theory why the energy levels become nearly equal. The repulsion of the 
nearly degenerate energy levels is explained by non-degenerate perturbation theory. We 
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show that the repulsion is caused by the opposite sign of the quadratic term in the 
perturbation series. This is confirmed by numerical calculations and is clearly visible in 
figure PB). 

In conclusion, we have shown that the present method can be used to calculate 
accurately the eigenvalues of the anharmonic oscillator in an external field. We also have 
shown that non-degenerate and degenerate perturbation theory are needed to explain 
the response of the energy of the ground state to an external field. 

In future work, we want to use these results to study the semi-classical behaviour 
of the quantum anharmonic oscillator. 
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